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Abstract 


In this thesis some weak external visibility problems have been considered. We have 
presented an optimal parallel algorithm for computing the shortest line segment from 
which a given convex polygon is weakly externally visible, and a sub-optimal parallel 
algorithm for computing the shortest line segment from which a given simple poly- 
gon is weakly externally visible. We also present an optimal parallel algorithm for 
computing the convex hull of a weakly externally visible simple polygon. The parallel 
computational model u.sed by our algorithms is the CREW PRAM. 
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Chapter 1 


Introduction 

1.1 Computational Geometry and Parallel Algo- 
rithms 

Computational geometry deals with geometric problems within the framework of de- 
sign and analysis of algorithms. It has applications in diverse areas such as computer- 
aided design, pattern recognition, computer graphics, VLSI design, image processing, 
vision, robotics, operation research, statistics, facilities planning and layout, etc. In 
the last two decades, computational geometry has attracted enormous research inter- 
est and is one of the fastest growing fields in computer science. This is because the 
knowledge and techniques stemming from research in this area has been successfully 
applied to many fields of science and engineering. 

Problems in the areas of artificial intelligence, medicine, energy resource manage- 
ment, metereology etc. require fast computing to process a vast amount of data, often 
in real-time. Over the past forty years, dramatic increanses in computing speed have 
been achieved. Most of these were due to the use of faster electronic components. 
Progress in this direction, however, is circumscribed by the laws of physics which put 
an upper bound on the speed of electronic devices. It appears that the only way 
around this problem is to use parallelism. The idea here is that if several operations 
are performed simultaneously, the time taken by any computation can be significantly 
reduced. Here, our computational tool is a parallel computer, that is a computer with 
many processing units or processors. Given a problem to be solved, it is divided into 
a number of subproblems. All of these subproblems are now solved simultaneously. 
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erich on a different processor. Tli<' results are then combined to produce an answer 
to tlie original problem. This is a very different model of computation from the one 
that has .served as the basis for the seciuenlial uniprocessor machines that have been 
built over the last forty years. 

Visibility problems are among the most fundamental ones in computational geom- 
etry. Visibility problems arise in many application areas, such as computer graphics, 
robotics, computer-aided design, facility layout etc. The focus of this thesis is on 
developing parallel algorithms for solving some weak external visibility problems. In 
particular, we are interested in computing the shortest line segment from which a 
polygon (simple or convex) is weakly externally visible. Initial results in this direc- 
tion were first obtained by Bhattacharya and Toussaint [9], who gave an optimal 
(sequential) algorithm for computing the shortest line segment from which a convex 
polygon is weakly externlly visible and then by Bhattacharya et al [8], who extended 
this result to the case of a simple polygon. 

1.2 Parallel Computational Model 

The commonly- accepted computational model for designing sequential algorithms is 
the Random .Access Machine(R.AM). The situation is somewhat different in parallel 
computing in that there is more than one accepted model of computation. However, 
one of the most commonly used ones is the Parallel Random Access Machine or 
PRAM, for short. The work described in this thesis is based on the PRAM model. 

PRAMs are also known as shared-memoiy'^ SIMD computers or array processors. 
An SIMD computer consists of some identical processors, all of which are supervised 
by a control unit. All the processors receive the same instruction broadcast from 
the control unit, but operate on different data sets from distinct data streams. For 
most problems, processors need to communicate among themselves during the com- 
putation in order to exchange data or intermediate results. This is done through a 
shared memory. Shared-memory SIMD computers can be further divided into three 
subclasses, depending on the policies used to resolve memory access conflicts. 
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1. Exclusive- Read, Exclusive- Write (EREW) SM SIMD computers : 

This model is the least powerful. In this, no two processors are allowed 
to read from or write into the same memory location simultaneousl 3 ^ 

2. Concurrent- Read, Exclusive- Write (CREW) SM SIMD comput- 
ers: In this model, multiple procesors are allowed to read from the same 
memory location but no two processors are allowed to write into the same 
memory location simultaneously. 

3. Concurrent- Read, Concurrent- Write (CRCW) SM SIMD comput- 
ers : This model is the most powerful. In this model, multiple processors 
are allowed to read and write from the same memory location simultane- 
ously. 

There are various methods for dealing with concurrent writes. These 
policies are, (i) all processors are allowed to write, provided that the value 
they are attempting to store are equal, otherwise access is denied to all 
processors, and (ii) the smallest numbered processor is allowed to write 
and access is denied to all other processors. 

In general, simulating a CREW PRAM or a CRCW PRAM algorithm on an 
EREW PRAM costs a slow-down b}'^ a logarithmic factor. Similarly, simulating a 
CRCW PRAM algorithm on a CREW PR.AM generally costs a slow-down by a log- 
arithmic factor. For more details on the PRAM model, see [1]. 

The PRAM model provides several advantages in the study of parallel algorithms. 
First, the PRAM seems to be suitable for studying the inherent parallelism of a prob- 
lem because in this model a researcher does not need to worry about interprocessor 
communication and processor synchronization, secondly, the PRAM can simulate any 
network model algorithm within asymptotically the same time and processor bound. 
Due to these advantages we use PRAM model to design our parallel algorithms. All 
algorithms in this thesis use the CREW PRAM model. 
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1.3 Our Thesis 

111 this thesis, we present new parallel techniques for solving some weak external 
visibility problems, related to simple and convex polygons. It is organized as follows. 
In chapter 2, we give a brief survey of some previous work, related to the results 
of our thesis. In chapter 3, we consider the weak external visibility problem on an 
n-vertex convex polygon C. 

We describe an optimal parallel cilgorithm for this problem which takes 0{log n) 
time, and uses njlog n processors. In chapter 4. we consider the same problem for 
a simple polygon. We give a sub-optimal parallel algorithm for this problem which 
takes 0[log^n) time, and uses 0{7iflog n) processors. Conclusion and future work is 
given in last chapter. 



Chapter 2 

Previous Work 


In this chapter, we give a brief survey of some previous work, related to the results 
of this thesis. 

Horn and Valentine [15] characterized L-sets in terms of their weak visibility prop- 
erties, and such characterizations for convex and stai’-shaped sets have been presented 
by Shermer and Toussaint [19]. 

Avis and Toussaint [5] showed that it can be determined in 0(n) time whether 
a given a simple polygon P is edge-visible from a specified edge. Sack and Suri [18] 
presented a linear-time algorithm for determining all the edges of P from which it is 
edge- visible. Yan Ke [16] considered the problem of detecting the weak visibility of 
a. simple polygon from an internal line segment, 1. He presented an 0{nlog n) time 
algorithm that tests if the polygon is weakly internally visible and reports one such 
segment when it is. He has also given an algorithm for the query version of this 
problem. Given a querj- line segment in P, whether P is weakly \'isible from it can 
be answered in 0{logn) time with 0{nlogn) preprocessing time aaid 0{n) space. He 
has also shown that a shortest such line segment can be found in 0{nlog n) time. 

Recently, Bhattacharya et al [7] considered the problem of determining the sector 
visibility of a polygon P. Informally, sector visibility addresses the question of exter- 
nal visibility along rays (or sight lines) whose angles are restricted to a sector(wedge) 
of specified angular width cr. Thus weak external visibility from a line corresponds to 
sector visibility, where a = w. They present algorithms which are 0(n) when <7 < tt, 
and <7 = TT but require n(?? log n) time when tt < a < 2 * tt. 
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A number of parallel algorithms have also been proposed for some other visibility 
problems, related to polygons. Atallah et al [4] presented an optimal parallel algo- 
rithm (on an EREW PRAM model) for the problem of computing the visible portion 
of a set of non-intersecting line segments from a point in the plane, using the ‘cas- 
cading divided-and-conquer’ technique. Another optimal parallel algorithm for this 
problem was gi\-en in [6]: this algorithm, however, assumes a CREW PRAM model 
and is based on the many- way divide- and-conquer technique. Chen [11] presented 
an optimal parallel algorithm for the case where the line segments ai’e the edges of a 
simple polygon. 

The kernal of a simple polygon P is the largest subset of P such that the en- 
tire polygon is visible from every point in that subset. Cole et al [12] presented an 
0{log n) time optimal parallel algorithm with 0{n/log n) CREW PRAM processors 
for computing this kernel. 

Goodrich et al [14] presented an 0{logn) time parallel algorithm with 0(n) CREW 
PRAM processors for computing the weak visibility region inside a simple polygon P 
from an edge of it. Goodrich et al [14] build a data, structure for ray shooting queries 
inside P in 0{log n) time using 0[n) CREW PRAM processors. The visibility graph 
on the vertices of P is obtained in 0{log n) time using O(nlog n -f k/log n) CREW 
PRAM processors [14], where k is the number of edges in the x'isibility graph. 

Chandra et al [10] presented a parallel algorithm for computing the complete 
\’isible polygon and the weak \'isible polygon inside a simple polygon P from a convex 
subpolygon C of P (a point p in P is completely visible from C if p is visible from 
every point of C, and is weakly visible form C if p is visible from at least one point 
of C); this algorithm takes O[logn) time and uses 0(n) CREW PRAM processors. 

Chen [11] presented an 0{logn) time optimal parallel CREW PRAM algorithm 
for computing all the edges of a given polygon P from which P is edge- visible, using 
0(n//opn) processors. 

Though the problem of weak internal visibility of a polygon has been well-studied, 
both in the sequential and parallel models of computation, hardly any result exists 
for the corresponding notion of weak external visibility of a polygon. Initial results 
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ill this direction wei-e first obtained by Bhattacharya et al [9], who presented a linear 
time sequential algorithm for computing the shortest line segment from which the 
given convex polygon is weakly externally visibile. Subsequently, Bhattacharya et 
al [8] extended this result to the case of a simple pol}^gon. 



Chapter 3 


Weak External 'Visibility of a 
Convex Polygon 

3.1 Introduction 

An 77-vertex coin-ex polygon C is wtukly txttrnally visible (wev, subsequently) from 
a line segment L, which lies outside C, if for every point .r on the boundary of C 
there is a point y on L such that the interior of does irot intersect the interior of 
C. In other words, an observer who traverses the line segment from one end to the 
other gets a complete view of the boundaiw of the polygon. Clearly, there could be 
many different line segments from which a polygon is weakly externally visible. For 
example, if the polygon is centrally symmetric there are an infinite number of such 
line segments. 

Bhattacharya and Toussaint [9] w'ere the first to consider the problem of weak 
external visibility of a convex polygon. They presented a sequential linear time algo- 
rithm for computing the shortest line segment from which such a pol3’’gon is weakly 
external^ visible. This was extended to the case of a simple polygon by Bhattacharya 
et al [8]. 

In this chapter, we consider the problem of designing an optimal parallel algorithm 
for the same problem in the CREW PRAM model of computation. 

The algorithm desci'ibed in [9] seems to be inherently sequential, but a couple 
of pertinent observations leads us to an optimal (parallel) algorithm that runs in 
0{log n) time and uses n/log n processors. 
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Figure 3.1: ray diagram 


The rest of this chapter is organised as follows. In section 2. we introduce some 
notations and definitions. Section 3 contains some results on a basic geometric min- 
imization problem and its extensions. Section 4 contains some useful observations. 
In section 5 we discuss a sub-optimal parallel algorithm, leading up to the optimal 
parallel algorithm, described in the last section. 

3.2 Preliminaries 

An n-veiiex convex polygon C is specified by a sequence (ci. C 2 , . . . , c„) of its vertices, 
in the order in which thej'’ are visited by a counter-clockwise walk along the boundary 
of C, beginning at ci. A chain on the boundary of C from c,-, say, counter clockwise 
to Cj,i ^ j is denoted by [ci^Cj]. An edge of C, joining Ci to Ci+i, is denoted by 
ti = CiCi +1 ( = c,+iC;), with the convention that c„+i = Ci. 

Translate each edge e,- so that the start vertex is at the origin of coordinates. 
What we get is the ray-diagram of the polygon, each sector representing one of its 
vertices. Further, the sectors occur in the same order as the vertices of the polygon. 
See Fig 3.1 

A pair of vertices of a convex polygon that admit parallel lines of support is called 
an antipodal pair. Given a vertex c,-, the set of vertices antipodal to it forms a chain 



10 


which is called the antipodal chain of c;. The initial and final vertices of this 
chain are found by drawing supporting lines to the polygon, parallel to the edges e,_i 
and fi respectively. 

The length of an antipodal chain is the number of vertices in it and is denoted by 

We denote a ray by op, where o is the origin of the ra}- and p an arbitrary point on 
it different from o; A cone is denoted by poq , where op and oq are rays with common 
origin o and the interior of the cone lies to right as we go from p to q via o, along the 
boundary of the cone. A line segment by pq, where p and q are the end points and a 
lint by where p and q are two arbitrary but distinct points on it. 

We denote by T’(C) the shortest line segment from which C is wev. 

3.3 Characterisation of the Shortest Line Seg- 
ment 

Our algorithm depends on a geometric minimization result and its variations proved 
in [9]. These are stated below without proof. 

Lemma 1: Let p be a point lying inside a cone W =^or and hk a line segment 
through p, with end-points h, k on oq and or respectively. Then \lik\ is a unimodal 
function of 0 tor (p < 9 < tt , where cp is the angle of the cone and 0 is the angle that 
hk makes with or. 

The above lenuna remains valid when we replace the point by a convex polygon 
and change other conditions in the following way. 

Lemma 2: Let C be a convex polygon contained in W, such that ^ and ^ 
ai'e tangent to it. Further, let hk be tangent to the subchain of C which is concave 
towards o. The end points h and k lie on oq and or respectively. Then j/ifc] is a 
unimodal function of ^ for <^ < ^ < r, where (p is the angle of the cone and 6 is the 
angle that hk makes with ^ 

Lemma 3: L*(C) has one of its endpoints on one of the bounding rays of 
and the other on the other bounding ray, for some value of i. 

Lemma 4: L*{C) is tangent to C. 
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Lemma 5: L’{C) is tangent to C at Cj and in c, 4 .ic,c,_i, then c, and Cj are 
antipodal to each other. 

Lemma 6: If x and y are two adjacent vertices then the last vertex in the 
antipodal chain of x(resp. y) is the first vertex of the antipodal chain of y{resp. x). 
Proof : It can be seen easilj' form ray diagram. Q.E.D. 

Thus for each vertex c,- there is an unique shortest line segment, which is tangent 
to its antipodal chain and has its end points on CtC,+i and So for ever}' vertex 

c, and for every vertex Cj in the antipodal chain of c,- we make use of lemma 2 to find 
the shortest line segment, which is tangent to the vertex cj and in c,+ic,c,_i. The 
shortest of all such line segments is the required line segment. 

3.4 Some Observations 

In this section we present some observations, useful to the development of our algo- 
rithm. 

Lemma 7: The total number of antipodal pairs in an n-vertex convex polygon is 
at most Zn/2. 

Proof : See Preparata and Shamos [17], Sec. 4.2.3. 

Lemma 8: The antipodal chain of any vertex can be computed in 0{logn) time 
using a single processor. 

proof : Since an antipodal chain is made up of a contiguous set of vertices, it is 
enough to locate the first and last ones, which can be done in 0[log n) (sequential) 
time by binary search. Q.E.D. 

3.5 A Sub-optimal Algorithm 

In this section, we present a basic procedure and a sub-optimal parallel algorithm. 

This basic procedure SHORT-LINE takes an antipodal pair (ci,Cj) as input and 
computes the shortest line segment, tangent to Cj, from which C is wev; the locus 
of the end-points of the line segment are determined by the edges of the polygon 
incident on c;. It repeats this by interchanging the roles ol c» and cj, and returns the 
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shorter of the two line segments. This procedure runs in 0(1) time by lemma 2. We 
describe it formally below. 


Procedure SHORT-LINE 

Input: An antipodal pair {ci,Cj) . 

The shortest line segment from which C is wev with respect 
to input antipodal pair. 


Output: 

step 1. 

Step 2. 

Step 3. 


Find the shortest line segment li , which is 


inside a cone c,+iC,Ci_i and tangent to C at Cj . 
Find the shortest line segment I 2 , which is 
inside a cone Cj+\cjCj^i and tangent to C at c,-. 
Return (min(/i, / 2 ) ) • 


It is easy to design a sub-optimal pai'allel algorithm, using the above procedure. 
Each vertex c,- is assigned a processor PEi, which computes its antipodal chain [cj, Ck]. 
By Lemma 8, this talces 0{log n) time. Depending on the length of this antipodal 
chain, we get three different cases. 

Case 1: |[Cj,Cjt]| = 1 i.e., j = k. 

Processor PE^ invokes the SHORT-LINE procedure with the antipodal 
pair {ci.Cj) as parameter. 

Case 2: \[cj,Ck]\ = 2 i.e., j + l = k. 

Processor PEi invokes the SHORT-LINE procedure twice, once with 
the antipodal pair (ci,Cj) as parameter and next with the pair (Ct,Cfc). 

Before discussing the third case we prove the following lemma. 

Lemma 9: Let the antipodal chain [cj,Ck] of a vertex c,- have length greater than 
2. Then the antipodal chain of any vertex, lying strictly between Cj and Ck, consists 
of Ci only. 
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Figure 3.2: 

Proof : We appeal to the ray-diagram for a neat proof of this. See Fig 3.2. 

To determine the antipodal chain of a vertex c,- from the ray-diagram, we extend the 
rays bounding its corresponding sector in this diagram in the opposite direction. This 
creates a sector. The vertices corresponding to the sectors that the sector so generated 
contain^ or intensect^constitute the antipodal chain of c,-. Therefore, the sector that 

we generate to determine the antipodal chain of a vertex that lies strictly between | ' 

I I 

the end-points of the chain, antipodal to c,-, is contained in the sector corresponding | i 

I ^ 

to c,-, which thus is its only antipodal vertex. This proves the lemma. Q.E.D. 

Now consider the case 3. 

Case 3: l[cj,c/ 5 ]| = m > 2 

By Lemma 9, the antipodal chain of a vertex in [cj+i, c^-i] consists 
of a only. Therefore, the allocation of work to the processors is done as 
follows. Processor PEi , where j < I < k , invokes the SHORT-LINE 
procedure on the antipodal pair (q, c,), while processor PEi taJces care of 
the antipodal pairs (c,-,Cj) and (q, c^-). 

Thus irrespective of the length of the antipodal chain of c,-, the associated processor 
invokes the SHORT-LINE procedure at most twice. 
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We give a formal description of this algorithm below. 

Algorithm SHORTEST- WEV-LIHE 
Input: An n-vertex convex polygon C . 

Output: The shortest line segment from which. C is wev 

step 1. Assign a processor to each vertex. 

Step 2. For i ;= 1 to n pardo 

a[i, 1] :=f irst vertex in the antipodal chain of Cj. 
a[i,2]:=last vertex in the antipodal chain of c; . 

Step 3. For i := 1 to n pardo 

h-SHORT-LINE(ci,Cj). 
h : =SHORT-LINE<icu c'k) • - 

l[i] :=minimum(/i, I 2 ) ■ 
endf or. 

Step 4. minimumd) . 

In this algorithm step 1 takes constant time. Step 2 takes 0{logn) time by lemma 
8. Step 3 takes constant time and step 4 takes 0{log n) time. Thus the shortest line 
segment can be found in 0{log n) time, using 0{n) processors. 

3.6 Optimal Parallel Algorithm 

In this section we present an optimal parallel algorithm, which takes 0{log n) time 
and uses 0{n./log n) processors. The essential idea of this algorithm is similar to 
that of the sub-optimal algorithm, presented in the previous section, with one major 
difference. Each processor is assigned to a group of log n contiguous vertices instead 
of one. 

We divide the n vertices into nflog n groups, Gi, G 2 , ■ • ■ iGn/iogm each consisting 
of logn contiguous vertices. Thus the group (?, contains the vertices from C(^i^iyogn+i 
to Ciiogn- Let processor PEi be assigned to group Gi. 
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Lach processor conij^utes the antij;odal chain of its group. That is, processor PEi 
finds the first vertex in the antipodal chain of and the last vertex in the 

antipodal chain of cn^gn- Let these be ajf?, 1] and ag[?',2] respectively. By lemma 6 
the antipodal chain of any vertex in Gi is a subchain of [aj,[z,l], 05(1,2]]. The length 
of the antipodal chain of a group can be 0(??), but we will show that it is sufficient 
to take care of 0(log n) antipodal pairs. 

Let a5[i,l] G Gj and ag[i,2] € Gk- Depending on the number of groups between 
Gj and Gjt, we get three cases as in the previous section. 


case 1 : j = ^M.e., ag[i,l] and 05[/‘,2] are in the same group. Then 
processor PEj invokes the SHORT-LINE procedure with antipodal pairs 
such as {p.q), where p G Gi and q G Gj. The following lemma provides a 
bound on the number of such pairs. 

Lemma 10 : Let «5[?',l],ag[i,2] G Gj. Then the number of antipodal pairs such 
^ {pi 9)1 where p G Gi and q G Gj is at most 2log n. 

Proof : We do an angular sweep of the ray-diagram in the counter-clockwise 
direction, with the initial position of the sweep line coinciding with the ray corre- 
sponding to the edge e^i-i)iogn+i- Each ray encountered in the sweep corresponds to 
a new antipodal pair, and since there are log n rays in the ray-diagram for any group 
the bound follows. Q.E.D. 

The time that this takes is obviously 0{log n). 


case 2: jf -t- 1 = A- i.e., 1] and ag[iy 2] are in two consecutive groups. 

Processor PEi invokes the SHORT-LINE procedure with antipodal pairs 
such as (p, q), where p G 6',- and q G Gj U Gk- The following lemma gives 
a bound on the number of such pairs. 

Lemma 11: Let 05(1,1] 6 Gj and a5[f,2] G Gk- Then the number of antipodal 
pairs (p, q) such that p G Gi and q E GjU Gk is at most Slog n. 
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Proof : Ihe proof is on the lines of the proof of lemma 10. We omit the details. 

Q.E.D. 


case 3: There is at least one group between Gj and Gk. 

It follows from lemma 9 that the antipodal chain of any group G in 
) ^. 7 + 2 ) • • • , Gk-i is a subchain of G',. Therefore the processor assigned 
to this group G invokes the SHORT-LINE procedure with its antipodal 
pairs, which by case 1 above is at most 2log n, so that it takes 0{log n) 
time. Processor P Ei invokes the SHORT- LINE procedure with antipodal 
pairs (p, q) such that p e G',- and q E Gj U Gk- By lemma 10 above, it 
takes Oilog n) time for this. 

Lemma 12 : Let ap[i,l] € Gj, a 5 [i, 2 ] € Gk- If there is at least one group strictly 
betM^een Gj and Gk, then the number of antipodal pairs such as {p,q) where p E Gi 
and q E Gj U Gk is at most Zlog n. 

Proof : The proof is on the lines of the proof of lemma 10. We omit the details. 

Q.E.D. 

Hence irrespective of the length of the antipodal chain of a group Gi, the associated 
processor invokes the SHORT-LINE procedure 0[logn) times. 

We now give a formal description of our optimal parallel algorithm below. 

Algorithm OPTIMAL- SHORTEST- WEV-Lim 
Input: An n-vertex convex polygon C . 

Output: The shortest line segment from, which C is wev 

step 1. Divide the n vertices into nflogn groups, 

each logn vertices. Let them call Gi,G 2 -, - . ■ iGn/iogn- 
Assign one processor to each group. 

For i := 1 to n/log n pardo 

«^[*,l]:=first vertex in the antipodal chain of . 

a^[?,2] :=last vertex in the antipodal chain of Cnogn- 


Step 2. 
Step 3. 
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Step 4. for i := 1 to n/log n pardo 

for each {p,q) such that p € G',- and q € Gj U Gk 
call SHORT-LINE(p,q) . 

endf or . 

Step 5. Find the minimum. 

Steps 1 and 2 in the above algorithm take constajit time. Step 3 takes 0{log n) 
by lemma S. Step 4 takes 0{log 7i) time by lemn:ias 10,11.12. The minimum can be 
computed in 0(log n) time. Hence our algorithm runs in O[log n) time. 



Chapter 4 

Weak External Visibility of a 
Simple Polygon 

4.1 Introduction 

111 this chapter, we discuss a. solution to the weak external visibility problem for a 
simple polygon. Bhattacharya et a/[8] were the first to consider the problem of weak 
external visibility of a simj^le polygon. They presented a (sequential) linear time 
algorithm for computing the shortest line segment from which a given simple polygon 
is wev. 

Their algorithm seems to be inherently sequential, but a couple of observations 
give us a sub-optimal parallel algorithm, which runs in 0{log^n) time using 0{nllogn) 
processors. The approach, however, is very different from that of the sequential 
algorithm. Henceforth by a polygon we shall always mean a simple polygon, unless 
otherwise stated. 

The chapter is divided into seven sections. In second section, we intz-oduce some 
notations, definitions and the rank tree data structure. The third section contains 
a characterisation of the shortest line segment from which a polygon is wev. In the 
fourth section we present a parallel algorithm for computing the convex hull of a 
polygon. The fifth section contains a discussion on envelope computation. In the 
sixth section we discuss how to compute a locally shortest line segment. In the final 
section we pi'esent our pai’allel algorithm. 
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4.2 Preliminaries 


111 this sect ion we present some notations, definitions and a data, structure called the 
runh iree, which is used lor efficient implementation of some of the parallel operations, 
required by our algorithm. 

4.2.1 Notation and Definitions 

An n-vertex simple polygon P is specified by a sequence (pi,p2,- • ■,Pn) of its vei'tices, 
in counter-clockwise order along the boundary. A portion of the boundary of a poly- 
gon P from Pi counter-clockwise to pj is denoted by [pi,pj]. .A,n edge of P, joining p,- 
to pi+i, is denoted by e,- = p,p,+i(— p,+ip,). By convention p„+i = pi. 

The convex hull of a polygon P is denoted by CH(P). The vertices of CH{P) 
are specified by the sequence (ci,C2,. . . ,c^-) of its vertices in counterclockwise order 
along the boundary. The angle made by a ray ^ with the positive x-axis is denoted 
by Q(op). 

A pocket of P is the region bounded by an edge CiCi+i and the subchain [c,',c,+i] 
of P. 

4.2.2 Rank Tree Data Structure 

The rank tree, data structure which supports the operations of parallel search, con- 
catenation and split, is described in[ll]; a similar data structure, called the hull tree, 
is described in[ 13 ]. 

A rank tree, RT{E), is a binary search tree with a set of edges E stored in its 
leaves in some specified order (e.g. by polar angle. If it stores points, then the points 
fire stored in order of increasing x-coordinates). Without loss of generality, assume 
that edges are stored in the leaves, and by increasing (or decreasing) order of polar 
angles. The leaves of RT{E) ai'e doubly-linked together. Each internal node v of 
RT{E) has four (^eH^ the first stores the number of leaves in RT(Ev), the second 
stores the edge t such that e has the smallest (or lai'gest) polar angle among the 
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edges stored in tlie leaves of the siid the other two respectively store the 

predecessoj- and the successor of e in the ordered set of edges. 

We denote the height of RT{E) (equal to the length of the longest root-to-leaf 
path in RT{E)) by h{RT(E)). Then in 0 {h{RT{E)) time. a processor can search in 
RT{E) for an edge with a given polai’ angle or search for the ranked edge stored 
in RT[E) (i.e., the leal of RT(E) in the left-to right order) using the first or the 
second feild stored in the internal nodes. For more details see[ll]. 

4.3 Characterisation of the Shortest Line Seg- 
ment 

Our algorithm depends on the generalisation of a certain geometric minimization 
result and some other related results proved in[S]. These are stated below without 
proof. 

Lemma 1: Let p be a point inside a concave cone W, formed by two polygonal 
chains oq = [oqiq2qz ■ ■ ■ 9m] arid or = [o7'ir2r3 . . . r„] and hk a line segment through p, 
with end-points /i, k on oq and or respectively. Then \hk\ is a unimodal function of 
the angle $ that hk makes with op; 9 € [a,/?], where a and j 3 are, respectively, the 
angles that the tangent s from p to or and op make with op. 

The above lemma remains \-alid when we replace the point inside the cone by a 
convex polygon and change othei' conditions in the following way. 

Lemma 2: Let C be a convex polygon lying inside a conca\-e cone W bounded by 
two convex polygonal chains oq — [09192 • • • 9m] and or = [o7‘ir2 . . . 7'm], with 0 being 
a vertex of C. If hk is a line segment lying inside the wedge, touching P, with its 
end points h and k lying on oq and 07' respectively, then \hk\ is a unimodal function 
of 9 , where 9 is the angle that hk makes with some ray op lying inside W; 9 € [a, 
where a and /9 are the angles that ^ makes with op when it is tangential to oq and 
or respectively. 

Lemma 3: The shortest line segment IF from which a polygon P is wev is tangent 
to its convex hull CH{P). 
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Let hk' be l.lie shoi't.est line segment from which a. polygon P is wcv. Assume that 
hk touches P at belonging to the antipodal chain of p. These two points split the 
boundary of P into two chains [p.z] and [r,p]. 

It is quite obvious that the vertices of the chain [p, *] are visible from hz only. 
This implies that all the^^external cones of support of the vertices belonging to the 
chain [p, z] intersect hz. Since hz is of minimum length, h lies on the upper envelope 
(with respect to z) of the CCW rays of the vertices of the chain [p, z]. We denote this 
upper envelope by CCW-envelope[p, z]. 

A similar argument shows that k lies on the upper envelope of the CW' rays of the 
vertices of the chain [z.p]. This upper envelope we denote by CW'-envelope[z,p]. 

In fact, one can make a stronger statement cis the following lemma shows. 

Lemma 4: Let [s,-,L] denote the antipodal chain of the vertex c,-, s and t be 
any tw'o vertices of [ 0 ,-,^^;] and [L-C;] respectively. The shortest line segment touching 
the convex chain [sj,L] with end points resting on CW-envelope[s,c,] and CCW- 
envelope[c,-, t] respectively is the shortest line segment antipodal to Ci from which the 
polygon is wev. 

The problem of computing the shortest line segment, therefore, is to determine 
the point of contact and the end points h and k. The point of contact is a vertex of 
the convex hull of P and the end points lie on a pair of CW and CCW envelopes as 
in the lemma above. 


4.4 Convex Hull Computation 


In this section we discuss a parallel algorithm for computing the convex hull of a 
wev polygon P. Let pi and p.,- be two vertices of P with the smallest and largest 
x-coordinate respectively. Clearly, pi and Pr are on CH[P). The line segment pip; 
divides CH{P) into two parts, an upper hull UH{P) and a lower hull LH{P). We 
will concentrate on the problem of computing UH{P), as the method for computing 
LH{P) is exactly similar. 

Given two disjoint upper hulls UH{R) and UH{S), the common tangent T, such 

of 

^ fc* r-, * * w 

U.T., 


that both UH{R) and^f/.^^^^re^b,4o^ft3&>s called the upper common tangent 

— t/A ■■ 
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U H{R) and UH(S). The following observations are useful in computing the upper 
hull. 

Lemma 5: Let pi and Pr be two vertices of P with the smallest and the largest 
x-coordinate respectively. A vertex p of the polygonal chain [pr,P/] that lies on or 
below pipr is not a vertex of UHiP). 

Proof : The proof follows from a convexity argument. Q.E.D. 

For any i,j and k such that r < i < j < k < I, assume that UH{\pi,pj]) and 
V H{\pj.pk]) have been computed and the common tangent T touches U H{\pi,pj]) at 
Pt^ and U H{[pj,pk\) at pt^. Since P is a wev polygon there exists a ray emanating 
from pj that does not intersect the polygon P and separates the polygonal chains 
[Pn,Pj] and 

Lemma 6: For any i,j and k such that r<i<j<k<l, assume that 
UH([pi,pj]) and UH{\pj,pk]) have been computed and given in rank trees RTi and 
RT 2 respectively. Then the upper convex hull UH{\pi,pk\) can be computed in time 
h = height{RTi) + height{RT 2 ), using a single processor. 

Proof : By using binary search, starting from pj oir the edges of upper hulls, 
we can find the upper common tangent T in time h = height{RTi) + htight{RT 2 )^ 
using a single processor. Assume that T touches I' pj]) at and U H{[pjiPk]) 
at ptj. Split the rank trees and merge the required rank trees: this can done in time 
proportional to h = height{RT\) + heighRRTi)- Q.E.D. 

Combining the abox’e observations with the optimal parallel algorithm for com- 
puting the convex hull of a sorted point set[13][ll], we can compute the convex hull of 
a wev polygon P in optimal time, that is, 0{d) time with 0{nld) processors. It must 
.be noted that we store the points corresponding to the vertices of our polygonal chain 
in the leaves of the rank tree in the order in which these appear along the boundary 
of the chain when we traverse it in, say, counter-clockwise order. 

4.5 Envelope Computation 

In this section we present an optimal parallel algorithm for computing envelopes. 
Before that, we present an algorithm for computing CW and CCW rays which is 
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required for envelope computation. 

4.5.1 Computation of Rays 

In this section, we present an optimal parallel algorithm for computing the external 
cone of support of each vertex of the polygon P. In other words, we compute the 
CW and CCW rays for each vertex. 

Chen[l 1] presented an 0{d) time optimal parallel algorithm for detecting the weak 
visibility of a. polygon from an edge using 0{n/d) processors, where d > log n. We 
make use of this algorithm. 

Using the parallel prefix technique, we divide the boundary of the polygon into 
groups of contiguous vertices such that : (i) the length of each group is greater than 
or equal to log n ; (ii) no pocket is divided in the middle i.e., the first and last vertices 
of any group are convex vertices. Assume that the group contains less than or 
equal to clog n vertices for some c < nflogn. We assign c processors to the group 
and compute the CW and CCW rays for each vertex in that group using Chen s 
algorithm. 

We formally sketch the algorithm below. 

Algorithm RAYS -COMPUTATION 
Input: An n-vertex simple polygon C . 

Output: CM' and CCW rays for each vertex p. 

Step 1. Assign 0(log n) contiguous convex vertices to each processor. 
Step 2. For i= 1 to k/log n pardo 

For j=(i-l)log n+1 to i log n do 

N[j]:=simple polygon chain length from cj to Cj+i 

steps. PREFIX-SUM (N). 

Step 4. For i= 1 to k/log n pardo 

For j=(i-l)log n+1 to ilog n do 

If ((trunc(N[j]/log n))<(trunc(N[j+l]/log n))) 
then MCj] := 1 else M[j] := 0. 
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steps. PREFIX-SUM (M) . 

Step 6. For i= 1 to k/log n pardo 

For j=(i-l)log n+1 to ilog n do 
If ^ M[j]) then 

Q[M[j],l] ;=trunc(N[j+l] /log n)-trunc(N[j]/log n) 
QCM[j].2] := j. 

Step 7. For i= 1 to k/log n pardo 

For j=(i-l)log n+1 to i log n do 

Assign Q[j,l] processors for the group 

of vertices from cq[j_i, 2 ] to 

Apply Chen's algorithm for each group. 


In the above algorithm step 1 takes 0(1) time while each of the steps from 2 to 
6 takes O(log n) time, using kllog n processors. Chen's algorithm takes 0{log n) 
time for each group and therefore the above algorithm takes 0{log n) time, using 
0{n/logn) processors. 

The rest of this chapter develops the computational machinery required for enve- 
lope computation. 

We divide the boundary of the polygon P into chains of length login) and assign 
a processor to each chain. Thus the polj^gonal chain [p^i-i^og n+iiP«/c>an] is assigned 
to the processor PEi. 

Let nii> 0) be the number of vertices of CH{P) in the polygonal chain which 
are indexed from zj to Processor PE^ finds the shortest line segment, which is 
tangent to the antipodal chain of convex vertices in its assigned chain. To do this 
each processor requires its own CW and CCW envelopes. We explain how to compute 
CCW-envelopes only, as the method for computing CW -envelopes is exactly similar-. 

The following observations come in handy. 

Let e-(= S7) be an edge of a CCW-envelope which is part of the ray CCW[z]. 
Such a ray is split into three parts, an initial par-t which precedes the edge e- and is 
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ccw[k] ccw[l] ccwfl] 



Figure 4.1: 


denoted by bef(ej)), the edge e’- and a final part, which is a ray, following the edge 
and denoted by aft(e'j). 

Lemma 7: Let the CCW-ray of vertex Pj+i, that is CCW[; + 1], intersect the 
CCW-envelope[l, j]. If CCW[t] is the last CCW-ray which contributed an edge to 
this envelope, where 1 < i < i, then a:(CCW[i]) > Q;(CCW[j -f 1]). 

Proof : Without loss of generality, we assume that the first edge in the CCW- 
envelope[l, j] is e[. From the definition of a CCW ray, all the vertices of polygon P 
lie to the left of CCW[1]. .Any vertex pk for k > i, is either to the left or to the right 
of the line(CCW[?]). 

If PA- is to the right of the line(CCW[i]) then CCW[A:] is bounded by the vertex p,. 
Therefore q(CX:W['/]) < q(CCW[A;]) and CCW[ii:] does not intersect CCW[i]. Thus 
CCW[A'] does not intersect the CCW-envelope[l,;]. See Fig 4.1(a). 

If Pa- is to the left of the line(CCW[?’]) then CCW[A:] intersects CCW[i] only if 
a{CCW[i]) > a(CCW[A:]). See Fig 4.1(b). 

Therefore for j > i, CCW[; 4-1] intersects the CCW-envelope[l, j] if Q(CCW[i]) > 
Q(CCW[i + Ij). Q.E.D. 

Froan the above lemma we can conclude that the set of CCW-rays of a polygonal 
chain which contribute an edge to the CCW-envelope of this chain are angularly 
ordered. A CCW-ray of a vertex in a polygonal chain is called effective if at any 
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liiiK' clurinjf IIk' iiu r<>ni('iital coiistrurtion ol tiio (.'('VV-eiivelope of this chain the ray 
coiitril)ut.ecl an edge to it. A subset oi the effective rays is a part of the final envelope. 
We can optimally compute effective rays by using the parallel prefix method. 

Lemma 8: Two CCW envelopes can intersect at most once. 

Proof : See[S]. 

Lemma 9: Let CCW-envelope[l,f] be kept in a rank tree RT. If CCW[j], for 
some j > z, is the next ray which intersects the CCW-envelope[l, i], then we can 
compute the CCW-envelope[l, j] in sequential-time 0(A), where h = heighi{RT). 

Proof : We use binary search to find the point of intersection between CCW[j] 
and CCW'-envelope[l, z]. Consider a ray CCW[A] such that k < z, such that is a 
edge of the CCW-envelope[l, z]. Let q be the point of intersection between CCW[A:] 
and CCW[j]. Depending on the location of q. we can eliminate a part of the CCW- 
envelope[l, z]. If q is on tf. then we have found the intersection between CCW[j] 
and CCW-envelope[l. /] (see Fig 4.2(a)). If q is on bef(e[,) then the intersection 
between CCW[j] and CCW-envelope[l,z] occurs on the part of the envelope before 
the edge e^. That is, we can eliminate the part of the envelope after the edge e). 
(see Fig 4.2(b)). Similarly, if q is on then the intersection between CCW[j] 

and CCV\^-envelope[l, z] occurs on the part of the envelope after the edge e'f.. That 
is, we can eliminate the part of the envelope before the edge e'/. (see Fig 4.2(c)). 
After finding the intersection, split the rank tree and merge with the ray CCW[_ 7 ]. 
The split operation can be done in time 0(A)[11]. Therefore we can compute the 
CCW-envelope[l,i] in 0(A) time with a single processor. Q.E.D. 

Lemma 10: Let (i?i, R 2 ) be two angularly ordered collection of CCW-rays. Given 
CCW-envelope(Ri) and CCW^-envelope(R 2 ) in rank trees RTj and RT 2 respectively, 
we can merge the two envelopes in sequential-time 0(/z). where h = height(RTi) -t- 
heightiRTi). 

Proof : By Lemma 8, two CCW''-envelopes can intersect at most once. So we 
can use double binary search, that is simultaneously use binary search on both the 
envelopes. We can compute the point of intersection in 0(/z) time. Therefore we can 
merge the two envelopes in 0{h) time with a single, processor. Q.E.O. 
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Lemma 11: If RT^ and RT 2 are rank trees corresponding to the envelopes CCW- 
( nv<;lope[l , ?'] and CCW-envelope[i + l,j], 1 < i < j < n, then we can merge the 
envelopes in time 0(/i), where h = htight{RT^) + heigkt{RT 2 ). 

Proof . The proof is exactly along the lines of the proof of lemma 10. We omit 
duplicating the details. Q.E.D. 

In the following lemmas we assume that we have more than one processor available 
to us. 

Lemma 12: Let us assume that the CCW-envelope[ij^i + l,?j] has been computed 
and given in Rlj lor j = I, . . . ,??? such that ij^i < ij and i„i < n with the convention 
that ?‘o = 0, then tor any j = I,...,??? we can construct a rank tree for CCW- 
euvelope[tj_i. /j]n CC'W-envelopefl.i,,,] in 0{h + logm) time using m processors, 

where h = max {height{RTj)} 

Proof : Our method for constructing the rank tree follows the treatment in[13]. 
Assign a processor to each pair (j, k){or k = 1, . . . ,y— 1, j+1 , ... ,m. Using lemma 11, 
we find the intersection between RTj and RTk . Let Ej^k be the edge on RTj containing 
the intersection point. It takes 0{h+logm) time [Oilogm) time is required to compute 
the value of h). Let Ej be the edge with smallest slope in {Ej^i.Ej^ 2 -. ■ ■ and 

Fj be the edge with lai’gest slope in {jEjj+i , Ejj +21 • • • j Fj and Fj can be found 

in 0(log m) time, using m processors. If the slope of Ej is greater then or equal to 
the slope of Fj then the part of the envelope from Ej to Fj is the required envelope, 
otherwise it is a null envelope. Q.E.D. 

Lemma 13: Let us assume that the CCW-envelope[ij_i+l, fj] has been computed 
and given in rank tree RTj for j = 1,. . . ,771 such that fj_i < ij and im < n with the 
convention that io = 0. Then the CCW-envelope[l,fm] can be computed using 
processors in 0{h + log 777.) time, 

where h = max {height{RTj)} 

l<j<m 

Proof : Assign m processors to each RTj and construct a rank tree RTj for 
J?r,nCCW-envelope[l, f„,] in 0{h + log m) time using lemma 12. Construct a rank 
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tree Rl', by building a. complete binary tree on top of RTj. It takes 0[log m) time. 
I'lierefore we can compute the CCW-envelopell,^,^] in 0{h + log m) time with 
processors. Q.E.D. 

Lemma 14: Let us assume that the CCW-enveIope[ij_i+l, ij] has been computed 
and given in RTj for j = 1, • • • , ni’ such that ij-i < ij and i-n < n with the convention 
that io = 0, then for j = we can compute the CCW-envelope[fj_i + l,im] 

using 0{m^) processors in 0{h + log m) time, 

whereh= max {height{RTj)} 

l<j<m 

Proof: Assign 0(rn.^} processors to each j = 1, .... m and construct the rank tree 
RTj for the CCW-envelope[fj_i,m] using lemma 13. Q.E.D. 

Our method of constructing CCW-envelopes is the following. We call procedure 
CCW-ENVELOPE. defined below, passing to it the set of CCW-rayes R and the 
integer log n, where n = |i2|. Our algorithm is bcised on lemma. 15 and the divide- 
and-conquer method of[ll]. 

Algorithm CCW-ENVELOPES 

Input: A set R of n CCW rays and an integer d. 

Output: n/log n CCW-envelopes. 

Step 1. If |J?1 < d then use one processor to compute the 
CCW-envelope . 

Step 2. If d < \R\ < d® then divide R into two subsets 
R^ and i?2 of equal size, and recursively solve 
the two subproblems in parallel. Then perform the merge 
operation on the output of the recursive call. 

Step 3. If > d® then partition R into g = {li?l/d)^^^ 
subsets Ri,...,Rg. Then recursively solve 
the g subproblems in parallel. Finally perform the 
computation using lemma 15, by using the output from 
the g recursive calls. 
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Lemma 15: The algorithm CCW-ENVELOPES takes 0{d) time using 0[nld) 
processors. It requiries 0[n?) space. 

Proof : See section 5.3.1 in [11]. 

4.6 Local Shortest Line Segment 

Let (c,, Cj) be an antipodal pair. Let denote its CCW-envelope and Ek its CW- 
envelope. 

To begin the search, an initial placement of the line segment hk is determined as 
followes: If a(c— TTcJ > • ISO - qIqZITq) then it can be determined by drawing a line 
through Cj parallel to and compute its intersection with Ek and Ek by using 

binary search method. Otherwise it is determined by the Cj_i, Cj. Similarly we can 
find the final placement of the line segment hk. 

The sequential search from an initial placement to the final placement is bounded 
by 0{n). That is, there may be at most 0{n) edges from an initial placement to the 
final placement on envelopes. But the total search time over all antipodal pairs is 
bounded by 0(n). If there are more than log n edges in the search range on the Eh 
then divide the envelope Ek into groups of logn and assign a processor to each gi'oup. 
Find the corresponding positions in Ek- VVe do the same for Ek- Each processor 
sequentially searches its respective range in Ek and Ek- 


Procedure LOCAL-SHORTEST-LINE 


Input: 

Output: 

step 1. 
Step 2. 


An antipodal pair {ci,Cj) and envelopes Eh and Ek- 
The shortest line segment from which P is wev with respect 
to input antipodal pair. 

Find the initial position and final position of hk . 
If number of edges on Eh between initial and final 
position is greater then logn then 

Divide the search range such that each 
subrange contains logn edges, assign 
a processor to each subrange and find the 
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corresponding position on Ek- 

Step 3. If number of edges on Ek between initial and final 
position is greater then log n then 

Divide the search range such that each 
subrange contains log n edges, assign 
a processor to each subrange and find the 
corresponding position on Ek- 
Step 4. Find the shortest line segment in each subrange 
sequentially. 

Step 5. Return the shortest line segment over all subranges. 

Lemma 16: Procedure LOCAL-SHORTEST-LINE takes 0{log n) time. 

4.7 Parallel Algorithm 

In this section we sketch the parallel algorithm. As in the case of the sequential 
algorithm, we divide our problem into two smaller subproblems. One subproblem 
is to find the shortest line segment for antipodal vertices of the vertices Cj, where 

j = 1,2, 771-1 .such that Ci is th antipodal vertex of c„j. The main idea is to 

divide a polygonal chain into subsets of O[logn) and assign a processor to each group 
and find the shortest line segment of wtv with respective to that group. 

Algorithm SHORTEST-LINE- SEGMENT 
Input: An n-veriex simple polygon P. 

Output: The shoHest line segment from which p is wev with 

step 1. Find the convex hull of given simple polygon CH{P). 

Step 2. Compute the CW and CCW rays for each vertex. 

Step 3. Compute the required nflogn CW and CCW- envelope s . 

Step 4. Divide the polygon chain into n/log n groups 

6\,G2,.--,Gnfic>gn and assign a processor to each group. 
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corresponding position on Ek- 

Step 3. If number of edges on Ek between initial and final 
position is greater then log n then 

Divide the search reoige such that each 
subrange contains log n edges, assign 
a processor to each subreinge and find the 
corresponding position on Ek- 

Step 4. Find the shortest line segment in each subrange 
sequentially. 

Step 5. Return the shortest line segment over all subranges. 


Lemma 16: Procedure LOCAL-SHORTE.ST-LINE takes 0{log n) time. 

4.7 Parallel Algorithm 

In this section we sketch the parallel algorithm. As in the case of the sequential 
algorithm, we divide our problem into two smaller subproblems. One subproblem 
is to find the shortest line segment for antipodal vertices of the vertices Cj, where 
j = 1,2, . . . .77? — 1 such that ci is th antipodal vertex of 0 ,^. The main idea is to 
divide a polygonal chain into subsets of 0{logn) and assign a processor to each group 
and find the shortest line segment of wev with respective to that group. 


Algorithm 
Input: 
Output: 
step 1. 
Step 2. 
Step 3. 
Step 4. 


SHORTEST-LINE-SEGMENT 
An n-vertex simple polygon P. 

The shortest line segment from which p is wev with 

Find the convex hull of given simple polygon CH{P). 

Compute the CW and CCW rays for each vertex. 

Compute the required n/logn CW and CCW-envelopes . 
Divide the polygon chain into njlogn groups 
Gi,G 2 ,.--,C?n/iosn and assign a processor to each group. 
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Step 5. For i:=l to n/log n pardo 

5.1. Find the antipodal chain of each group. 

Let Ai be the antipodal chain of group Gi 
(i.e., first antipodal vertex of first 

convex vertex aoid last antipodal vertex of the 
last convex vertex in that group.) 

5.2. If the aintipodal chain length of a group is less than 
log n then find the shortest line segment sequentially. 

5.3. If the aintipodal chain length of a group is greater than 
log n then 

5.3.1. Divide the polygon chain Aj into groups 

each logn length and assign a 

processor to each group. 

5.3.2. Find the antipodal chain of each polygon chain 

Aij . Let them call G',j , Gi^ 

5.3.3. Compute the CW and CCW envelopes for each group 

5.3.4. compute the shortest line segment for each group. 
Step 6. Find the shortest line segment. 


Each processor is assigned to two polygonal chains G,-^ and Ai^, which are an- 
tipodal to each other. The length of these chains are at most log n. Each processor 
computes the shortest line segment, which is tangent to the convex polygon CH{P) 
at some vertex in the polygon chain .4,-^. So the maximum number of processors 
required is 0{n/log n) 

In the CCW-ENVELOPES algorithm we divided the polygonal chain into groups 
of log n. In each group the first vertex need not be a convex vertex. We need to 
adjust these envelopes so that the first vertex is convex. We make this adjustment, 
using the merge and unmerge operations. We need to perform merge operation on the 
CW-envelope and unmerge operation on CCW-envelope. By Lemma 11 a merging 
operation can be done in 0{log n) time using a single processor. For an unmerging 
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operation, let us assume that jirocessor PE is assigned to tlie polygonal chain G = 
and it is required to unmerge the C('W-envelope[/, j], where i < j < k, from 
its CCW-envelope[?., ??.]. Processor PE sequentially computes the CCW-envelope[j', A:] 
and merges it with the CX^W-envelope of the next processor. The length of the 
polygon chain G is less than log ii. So the CC\'V-envelope[_/, k] can be constructed in 
0{logn) time. Therefore, unmerging can be done in 0{logn) using a single processor. 

In the above algorithm each of the steps from 1 to 3 takes 0{log n) time using 
0{nllogn) processors. We can compute the antipodal chain of any contiguous group 
in 0{log n) time. So steps 5.1 and 5.3.2 can be performed in 0{log n) time. We can 
compute the envelopes for the new groups in 0{log n) time using merge operation. 
So step 5.3.3 can be computed in 0{logn) time. 

Each processor should find the shortest line segment with respective to the as- 
signed groups. Each processor takes one antipodal pair at a time from its groups and 
calls the LOCAL-SHORTEST-LINE In the LOCAL-SHORTEST-LINE 

procedure we have divided the search range and assigned the processors to each range. 
The total search range is at most 0{n). Therfore the maximum number of processors 
required is 0{nflog n). So we can perform steps 5.2 and 5.3.4 in 0{log n) time. 

We need to perform merge and unmerge operation for each convex vertex in its 
assigned groups. Each group has a maximum log n convex vertices. So it takes 
0{log‘^n) time. 

Lemma 17: Algorithm SHORTEST-LINE SEGMENT computes the shortest line 
segment from which a given simple polygon is wtv in 0{log^n) time using 0{njlogn) 
processors. 



Chapter 5 
Conclusion 


In this thesis, we have presented parallel algorithms for solving some weak exter- 
nal visibility problems, related to simple and convex polygons on CREW PRAM. 
The basis of these algorithms are very different from those used for designing the 
corresponding sequential algorithms. 

I 

In chapter 3 we have designed an optimal parallel algorithm for computing shortest 
line segment from which given convex polygon C is weakly externally visible. The 
algorithm is based on some properties of antipodal pairs of points. These are: (i) If 
the antipodal chain length of a vertex is greater then two, then all vertices in that 
chain, except for the first and the last, have an antipodal chain of length exactly 
one. (ii) The number of antipodal pairs is 0(7?). Using these observations, we have 
designed an 0{logn) time parallel algorithm that requix'es 0[nllogn) CREW PRAM 
processors. 

We have presented an optimal parallel algorithm for computing the external cone 
of support of a vertex p of a given simple polygon P. The algorithm returns this cone 
if P is weaklj^ externally visible else retuns that P is not weakly externally visible. 
We use a techinique presented in [11]. Our algorithm takes 0{log n) time, using 
0{nflogn) CREW PRAM processors. 

We have presented an optimal parallel algorithm for computing the convex hull 
of given weakly externally visible simple polygon P. The observation involved in our 
algorithm is that for every vertex p on the boundary of P, there exists a ray which 
emanates from p and does not intersect the boundary of P. Using this observation, 
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we have reduced our problem to the pi oblem of computing the convex hull of a sorted 
set of points [11] [13]. 

We have also presented a sub-optimal parallel algorithm for computing the shortest 
line segment from which given simple pol)'gon P is weakly externally visible. We 
use the rank tree data structure of [11] and a divide-and-conquer technique. Our 
algorithm takes 0{log^n) time, using 0{n/log n) CREW PRAM processors. 

This research presented in this thesis can be extended in the following directions. 

One direction is to look into the design of an optimal parallel algorithm for the 
problem, discussed in the previous chapter. This essentially boils down to the problem 
of merging and unmerging the envelopes in linear time. 

Another direction that can be explored is to see whether the approach presented 
here can be used to design an optimal parallel algorithm to compute the shortest 
internal line segment from which a polygon is weakly internallj'^ visible. 

A third direction is to study these problems in higher dimensions. Designing fast 
sequential algorithms for these problems in higher dimensions are still open. 
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